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Motivated by the observation of several molecule candidates in the heavy quark sector, we discuss 
the possibihty of a state with J^'~^ = S""*". In a one-boson-exchange model investigation for the S 
wave C = + D*D2 states, one finds that the strongest attraction is in the case J = 3 and 7 = for 
both TT and a exchanges. Numerical analysis indicates that this hadronic bound state may exist. If 
a state around the D* D2 threshold (Ri4472 MeV) in the channel J /%l)uj (P wave) is observed, the 
heavy quark spin symmetry implies that it is not a cc meson and the J^^ are very likely to be 3 ''. 



PACS numbers: 12.39. Pn, 12.40.Yx, 13.75.Lb, 14.40.Rt 



I. INTRODUCTION 

Mesons with exotic properties play an important role in understanding the nature of strong interactions. The 
observation of the so called XYZ states in the heavy quark sector has triggered lots of discussions on their quark 
structures, decays, and formation mechanisms. It also motivates people to study new states beyond the quark model 
assignments. 

The X(3872), first observed in the Jipn^n^ invariant mass distribution by Belle collaboration in 2003 [l[, is the 
strangest heavy quark state. Even now, its angular momentum is not determined. Since its extreme closeness to 
the D°Z)*° threshold, lots of discussions about its properties are based on the molecule assumption. However, it is 
very difficult to identify the X(3872) as a shallow bound state of D°Z)*° since there arc no explicitly exotic molecule 
properties. 

A charged charmonium- or bottomonium-like meson labeled as Z is absolutely exotic because its number of quarks 
and antiquarks must be four or more. Such states include the .Z(4430) observed in the ip'n^ mass distribution [2], the 
Zi(4050) and ^2(4250) observed in the Xci'^~^ mass distribution [3, and the ^6(10610) and ^^(lOeSO) in the mass 
spectra of the T(niS')7r^ (71=1, 2, 3) and hi,{mP) (m=l,2) They are all observed by Belle Collaboration. Though 
BABAR has not confirmed them 0|, the existence signal of multiquark states is still exciting. Since Z(4430) is 
around the D*Di threshold, Zh(10610) is around the BB* threshold, and Zb(10650) is around the B*B* threshold, 
molecular models seem to be applicable to their structure investigations [7Hl3l| . 

To identify a state as a molecule is an important issue in hadron studies. One should consider not only bound state 
problem of two hadrons, but also how to observe a molecular state in possible production processes. In Refs. [l3-[l3i 
bound states of and "ScD* were studied. Since the quantum numbers are the same as the nuclcon but the masses 
are much higher, identifying them as multiquark baryons is rather apparent. To obtain a deeper understanding of the 
strong interaction, it is necessary to explore possible molecules with explicitly exotic quantum numbers. 

Quark model gives us a constraint on the quantum numbers of a meson, namely, a meson with J^^ ~ , 0"' , 
1 '', 2"' , 3 '", • • • could not be a qq state, but it may be a multiquark state. So the study on such states may deepen 
our understanding of nature. If two qq mesons can form a molecule with such quantum numbers, one gets the simplest 
configuration. Next simpler configuration is the baryon-antibaryon case. A possible place to search for them is around 
hadron-hadron thresholds. There are some discussions on low spin heavy quark exotic states in Refs. [l^, [l^. Here 
we would like to discuss the possibility of a higher spin state, J^^ = 3 ^ . One will see that identification of it from 
strong decay is possible. 

First, we check meson- antimcson systems that can form 3 states, where meson (antimeson) means that its quark 
structure is cq (cq). The established mesons may be found in the Particle Data Book [2^. One checks various 
combinations and finds that the lowest S-wave system is D*D2- The next S-wave one is D*D*2- Between these two 
thresholds, one needs D or G wave to combine other meson-antimcson pairs (see Fig. [l}. Below the threshold of 
D*D2, the orbital angular momentum is D, F, or G-wave. Above the £'*D*2 threshold, a partial wave of P, F, or 
H is needed. Since the difference between these two thresholds is more than 200 MeV, one may neglect the channel 
couphng and choose the D*D2 system to study. 
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FIG. 1: Thresholds of J^^ = 3 meson-antimeson systems between that of D* D2 and that of DlDl2- S, D, G, and I are 
orbital angular momenta. 



Secondly, we check baryon-antibaryon systems. If one combines the established cqq baryons and their antibaryons, 
one finds that the lowest S-wavc threshold is for Ac(2880) and Ac (w 5168 McV). Even for the lowest threshold of 
Ac(2595) and Ac in F-wave, the value (f* 4879 MeV) is still higher than that of D*D*2. Thus, we may safely ignore 
the possible baryon-antibaryon contributions in this study. 

In a 3 D*D2 state, partial waves of S, D, G, and / may all contribute. As a first step exploration, we consider 
only the dominant S-wave interactions. Possible coupled channel effects will be deferred to future works. The present 
study is organized as follows. After the introduction in Sec. U we present the main ingredients for our study in Sec. 
im Then we give the numerical results in Sec. IIIIl The final part is for discussions and conclusions. 



II. WAVEFUNCTIONS, AMPLITUDES, AND LAGRANGIAN 



We study the meson-antimeson bound state problem in a meson exchange model. The potential is derived from 
the scattering amplitudes [2l| and the fiavor wave functions of the system are necessary. Since the states we are 
discussing have a definite C-parity while the combination of cq and cq mesons does not, a relative sign problem arises 
between the two parts of a flavor wave function. One has to find the relation between the fiavor wave function and 
the potential with definite C-parity. There are some discussions about this problem in the literatures |2ll425l. Here 
we revisit it by using the G-parity transformation rule which relates the amplitudes between NN and NN |26| . 

Since D mesons do not have defined C-parity, we may assume arbitrary complex phases a and /? under the C-parity 
transformations 
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According to the SU(2) transformation, one finds the following isospin doublets 
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from which the G-parity transformations read 
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Similar to the study of the D*Di bound state problem l2l|, one may construct several states from D* and 1)2 ■ 
Here we concentrate only on the C — + case. If the system is an isovector (isoscalar), we label it Zj {Xj) where J is 
the angular momentum. Explicitly, one has the G-parity eigenstates 
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+cai/3l{D;- D*+ + D*+ D;-) + ca2Pl{D;°D*" + D*"D;°) , (4) 

where c = 1 is the C-parity and the superscript indicates the electric charge. One may check 

GZ° = -cZ'j, CX° = cX^j. (5) 

The procedure to derive the potential is similar to that in Ref. [2l|. Now we calculate the amplitude T{Zj) = 
{Zj\T\Zj) with the G-parity transformation rule (|3]). We just present several terms to illustrate the derivation. 
Together with the above wave function, one has 

= — |7'[D*+^D*+.r>*"^r>*"] + "i"2/3i/32^[D*+^D*M)-"^r>*+] + c^i^'2P^P2T[di+ ^d'+ ,D'»^d'^°] ^ | (6) 

In fact, the convention aia\ = /3i/3| is implied in the Lagrangian in Eq. So aia\f3\l32 = a\a2l3i(i\ ~ 1 and one 
finally gets 

"^(^j) = -jG"" {^[L>' + -^D* + ,_D'0->r>*0] ^ ^[Z5* + ^_DJ'\D*"^_D' + ] ~ ^[DjO-^_D*+ ,_D' + ^D*0] + ^[-Dj"^r'3*",D' + ^D* + ] 

It is obvious that we may calculate the potential of meson-antimeson interaction from that of meson-meson together 
with a given Lagrangian for (oq) meson fields. The arbitrary relative phase in the flavor wave function of a meson- 
antimeson system is canceled in this procedure. To derive the explicit expression of the potential, one needs interaction 
Lagrangian. 

The Lagrangian for pion interactions in the heavy quark limit and chiral limit reads [27l [28j 
= .gTr[i?475i7] + .g"Tr[r^475r^] 

+{^Tr[T^(I?^4)75i/] + /i.e.} + {^Tr[r^(#l^)75i?] + /i.e.}, (8) 

where 

a = ^[^*^7m + ^75], ff = 7"i?V 

The flelds P* = !?*+), and = (D'^ ,D*^) annihilate the eg mesons. P and P\ have similar form but we do 

not involve them in the following calculation. The axial vector field is defined as ~ ^{^d^^ ~ ^d^^) with 
C = exp(i7W//), / = 132 MeV and 



-71/ 

If we further consider a exchange, one needs additional interaction terms 

= .gaTr[i/<7i/]+.g;'Tr[T'^af'^] + :^Tr[r^(a^a)i7 + i?(9^(T)T^]. (11) 

The coupling constants must be determined in order for numerical analysis. One extracts the pion coupling constant 
g from the decay D* Dtt: g = 0.59 ± 0.07 ± 0.01 [H. For = -^^J^, we use the value 0.55 GeV"^ estimated in 
Ref. pTf . To determine the coupling constant g", we turn to the chiral quark model (30| with which one may get the 
relation g" = —g. 
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For the sigma coupling constants, we can just get estimates from the chiral quark model or the chiral multiplet 



assumption [31[. These approaches have been used in the baryon case |32| for the purpose of cross checking, where 
we get consistent results. Now the former method may give the relation g'^ = —g^ and the value g^ ~ gch — 2.621 if 
one adopts the Lagrangian [s^ 

Li = -,9c/i'0(o" + i757raT-a)'0, (12) 

where tp = (u, d)^ is the quark field and Ta the Pauli matrix. One should note the normalization problem in this 
approach [28l. |33|. However, if one estimates g^ from the chiral multiple assumption, a value less than 1 is obtained. 
For the remaining h'^, no available approach may be used. Since the large uncertainties of the coupling constants, we 
will select several values to see the u-exchange effects on the conclusions. 

In deriving the above relations for the coupling constants, we have used the polarization vectors e'^^ = -^(0, ±1, i, 0) 

and Eg = (0, 0, 0, —1) for the vector meson D* and 



= A/i[4i£-i+£-i£';i+2e[;£S]. (13) 



for the tensor meson Z?2 [34 



III. POTENTIALS AND NUMERICAL ANALYSIS 

Now one may derive the potentials through the amplitudes in ([7]). Using the same procedure as Ref. [2l|, one gets 
the one-pion-exchange potential (OPEP) for S-wave interaction in the case 1 = 1 



15/2 
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(14) 



where /i = {rriD^ — mjj* — , and 



-1, (J = 3) 



i (J = 3) 



a,. = i i, (J = 2) , Ce = { -|, (J = 2) 

I, (^ = 1) I -I, (^ = 1) 



(15) 



There are two parts in the potential: direct part and spin-exchange part. The later corresponds to the terms containing 
c in Eq. 0. For the case of / = 0, ViXj) = -W{Zj). 
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The singular behavior at small distances needs to be regularized [35| . If a form factor FF 
each vertex, one finally has 
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where r] = A/A2"^^7ri2^, and a = h? — (muj — mu*)^. 
Similarly, the one-cr-exchange potential (OsEP) is 
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FIG. 2: One-pion-exchange potentials for (a) Z states and (b) X states with the cutoff A = 1 GeV. 
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FIG. 3: One-sigma-exchange potentials for X and Z states with the cutoff A = 1 GeV and the coupling constants g'J. = 
-1.0, h' = 1.0. 
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2, and ^ for J 
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1. 



where fia = — (m_D2 — "^d*)^, '7cr = a/A^ — . The coefficient = 1 for J = 3, ^ for J = 

The spin-dependent nature of OsEP conies from the third couphng term in the Lagrangian (jlip. 

Before the numerical calculation, we take a look at the relative strengthes of potentials. For the meson masses, we 
use = 137.27 MeV, mo* ^ 2008.63 MeV, and mo^ = 2463.5 MeV 0. We plot OPEPs with A = 1 GeV in Fig. 
[2j It is obvious that X3 is the most attractive case. For the cr meson exchange contribution, we use ma- — 600 MeV 
for the illustration. In Fig. [3l we show OsEPs with g'l = —cja = —1.0, h'^ = 1.0, and A = 1 GeV. It is interesting 
that the potential for is also the most attractive one. Thus the long-range and medium-range meson-exchanges 
are both helpful for the formation of a l'~^{J^'~^) = 0^(3 ^) state. 

Now we turn to the numerical results for the OPEP case by solving the Schrodinger equation. In the potential, there 
is an unknown phcnomcnological cutoff parameter A. It incorporates the size information of the interacting mesons. 
If A goes to infinity, the potential describes the interactions of structureless mesons. A small cutoff is relevant to the 
real case. In principle, an appropriate value should be around 1 GeV which is realized from the nuclear potential 
models. Since we do not want to give accurate prediction on the binding energy, this parameter is not fixed. We 
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just tune its value in a range and check whether a bound state exists or not. The resuhs of binding energy (B.E.) 
and root-mean-square radius (r^ms) for with various A are presented in Tab. HI where the cases for r^ms <0.8 
fm or A >4 GeV are neglected. Similarly, one may get numerical results for other possibilities, which are also given 
in that table. The resultant cutoff much larger than 1 GeV indicates that the attraction is not strong enough for 
the formation of a hadronic bound state. Comparing the three cases in the table, of course X3 is more likely to be 
existent. 



State A (GeV) B.E. (MeV) r™, (fm) 



V3 


2.3 


0.6 


3.8 




2.4 


3.7 


1.6 




2.5 


9.9 


1.0 




2.6 


19.8 


0.8 


Z2 


3.8 


2.5 


1.9 


Zi 


3.6 


1.9 


2.2 




3.7 


8.4 


1.1 



TABLE I: Cutoff (A), binding energy (B.E.) and root-mean-square radius (r^ms) for X and Z states with OPEP. We do not 
show results for A > 4 GeV or Vrms < 0.8 fm. 

The minimal cutoff for a binding solution is a little larger than 2 GeV if we consider only 7r-exchange. This means 
that additional attraction may lower the value to a more appropriate number. We would like to check how much 
attraction the sigma meson contributes. Because of the large uncertainty of the coupling constants, we take three 
sets of them: (1) go.=2.621, g'^ = —go, h'^ = 0, which corresponds to neglecting spin-exchange potential; (2) g^y = 1.0, 
= —.go-, ~ 1; ^-i^d (3) ga = 2.621, g" = —g^, h'^ = 2.621. The last set is the most attractive case. After the 
solution of the Schrodingcr equation, the cutoff parameters satisfying the condition r^ms >0.8 fm and A <4 GeV are 
summarized in Tab. |TI1 From the resultant cutoff parameters, we see that the existence of X3 is probable. Its mass 
should be around the D* threshold (w 4472 MeV). 



States 


Set 1 


Set 2 


Set 3 


V3 


1.7-2.2 


1.5-1.7 


1.0-1.1 


X2 






1.4-1.5 


Xi 








Z3 




2.7-3.0 


1.1-1.2 


Z2 


2.6-3.2 


2.5-2.7 


1.2-1.3 


Zi 


2.2-2.9 


2.8-3.1 


1.5-1.7 



TABLE II: Cutoff values (GeV) for X and Z states with OPEP+OsEP when binding solutions exist. We do not show cutoffs 
if A > 4 GeV or rrms < 0.8 fm. 



IV. DISCUSSIONS AND CONCLUSIONS 

From the meson exchange potentials and the numerical analysis, one has found that the most probable molecule 
in the C = + D*D2 system is X3. If the state really exists, it may decay through its components, i.e. D* — > Dn, 
D2 Dtt, or D2 — >■ D*Tr. The X3 may also decay through the quark rearrangement, i.e. the final states are a cc 
meson and a qq {q ~ u, d) meson. The later type decay may be used to identify the exotic quantum numbers. Here 
we focus only on this case. 

For convenience of discussion, we assume that L is the relative orbital momentum between the cc and the qq mesons 
and relax the isospin requirement temporarily. Since the spins of the charm quark and the light quark in both D* 
and D2 are parallel, the spin of cc in V3 must be 1. According to the heavy quark spin symmetry, the spin of the 
final charmonium after rearrangement should also he S = 1. Thus the final cc state can only be ip or Xcj- The decay 
channels are obtained as follows: 

(1) If the final cc is J/tjj, the J^'^ of the produced qq meson may be (1 — 5)~~ for L = 1, (1,3,5)"' for L = 2, 
(1 — 7) for L = 3, and so on. After some inspections on the meson masses, one finds that kinematically allowed 
decays for the X3 molecule are just J/ipp and J/^puj with L ~ 1,3, 5, and J/tphi{1170) and J/V'&i(1235) with L = 2, 4. 

If it is ip{2S), the kinematically allowed decays are ip{2S)p and ip{2S)uj with L = 1, 3, 5. 
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(2) If the cc is Xco, the J^^ of the qq meson may be (2 - 4)++ for L = 1, (2, 4)"+ for L = 2, (0 - 6)++ for L = 3, 
and so on. The kinematicahy allowed decays are Xco/o(500), Xco/o(980), and Xcoao(980) with L = 3. 

(3) If the cc is Xci, the J^^ of the qq meson may be (2,4)"+ for L = 0, (1 - 5)++ for L = 1, (0,2,4,6)"+ for 
i = 2, (0 ~ 7)++ for L = 3, and so on. The kinematically allowed decays are Xci^r, XciVi XcW with L = 2,4, and 
Xci/o(500) with L = 3. 

(4) If the cc is Xc2, the J^'^' of the meson may be (2,4)-+ for L = 0, (0 - 6)++ for i = 1, (0,2,4,6)"+ for 
L = 2, (0 ^ 8)++ for L — 3, and so on. The kinematically allowed decays are Xc2/o(500) with L = 1,3,5, and Xc2'^ 
and Xc2?? with L = 2, 4. 

Therefore, the allowed two-body strong decays for X3 arc J/tpuj (PFH). -0(25')^ (PFH), J/'(/;/ii(1170) (DG), 
Xco/o(500) (F), Xco/o(980) (F), xaV (DG), XcW (DG), Xci/o(500) (F), Xc2/o(500) (PFH), and Xc2V (DG). There is 
no S-wave decay. Because high L processes are suppressed and ■0(25') and Xc2 are excited states, the simplest way to 
identify may be through the J/ipuj channel. 

Let us analyze the J^^ of an assumed state X(4472) observed in the J/ipur mass distribution. Since J/ijj and u 
are both J^'~' = 1"" mesons, the quantum numbers of J/ipu! are (0, 1,2)++ for S-wave combination, (0 ~ 3) ^ for 
P-wave combination, (0 ~ 4)++ for D-wave combination, and so on. What we are interested in is the case that the 
partial wave is determined to be P. If X were a conventional cc meson, the state is 77c (4472) and the spin of cc must 
be 0. Because of the heavy quark spin symmetry, the decay 77c (4472) J/ipuj is suppressed. Then X(4472) could be 
a hadronic state. Although other meson-antimeson pairs may also form molecules with J^'~' = (0 ~ 2) ^, the masses 
are smaller. Therefore, based on our numerical analysis, this X around the D*D2 threshold is very likely to be a 
state with the exotic J^'~^ = 3 ^. 

If one wants to look for Z3, one can use those kinematically allowed decay channels, J /ipP (PFH): i^{'^S)p (PFH), 
J/0&i(1235) (DG), Xcoao(980) (F), Xciti" (DG), and Xc27i" (DG). The practical way to identify the J^^ is to analyze 
the partial wave of J /"ipP- The search is also helpful to test the meson exchange models. 

We consider only S-wave interactions of D* and D2 in this paper. Higher partial waves also have contributions to 
the 3 ^ molecular bound state. Additional attraction from the channel coupling may reinforce the present preliminary 
result. Future studies on such effects and production cross section may be helpful for the search on exotic states. 

Replacing a c quark with a h quark, one may study the bottom case. Because the production of a hidden bottom 
molecule B* B2 needs much higher energy and the production cross section is smaller, it is difBcult for experimentalists 
to explore this case in near future. However, with the replacement c — >■ s, one may study whether there is a bound 
state or resonance with J^^' = 3 ^ near the K*K2{~ 2322 McV) threshold. If such a state exists, it may decay into 
ujcf) and can be detected. 

In short summary, we have investigated whether hadronic bound states exist in the D* D2 system in a one-boson- 
exchange model. The C = + case is discussed in this paper. We find that the I^{JP^) = 0+(3-+) X3 state is the 
most probable one. A feasible place to identify it may be in the invariant mass distribution of J/ipu around 4472 
MeV. A similar study for a state around 2322 MeV in the cj0 mass distribution is also called for. 
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